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The observation of 1.97 ± 0.04 solar-mass neutron-like star gives constraint on the equation of 
state (EOS) of cold, condensed matter. In this paper, the EOS for both pure quark star and hybrid 
star with a quark core described by quasi-particle model are considered. The parameters of quasi- 
particle model which affect the mass of both quark star and hybrid star can be constrained by the 
observation. 
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I. INTRODUCTION 

It is well known that the equation of state (EOS) for 
cold and dense strongly interacting matter plays a key 
role in the study of the neutron stars in astrophysics 
[l], 0. The neutron star are presumed to contain the 
densest matter in the cosmos, which provide a natural 
laboratory for cold, condensed matter. The recent obser- 
vation [I of the pulsar PSR J1614-2230 with a mass of 
1.97±O.O4M0 (Mq denotes the mass of the sun) gives 

a strong constraint on the EOS of strongly interacting 
matter. 

Just as the discovery of neutron led to the idea of 
neutron star, the conception of quark naturally stimu- 
lated many physicists and astrophysicists to suggest that 
neutron-like star may be quark or hybrid star (Quark or 
hybrid star is composed, in whole or in part, of quark 
matter) At present, it is still not possible to get 

a reliable EOS of quark matter from the first princi- 
ple of Quantum Chromodynamics (QCD). In this case, 
a commonly used approach in the study of high den- 
sity quark matter is to give up on first principle cal- 
culations and resort to phenomenological QCD models, 
such as MIT bag model EMI, NJL model [TUH, P er " 
turbative QCD model [19|-|22| and quasi-particle model 
[23U2II ] . Among these models, the quasi-particle quark- 
gluon plasma model with a few fitting parameters, has 
been widely used to reproduce the properties of the QCD 
plasma at finite temperature and zero chemical potential 
[I2I [27-33]. However, due to lack of experimental data 
on dense strongly interacting matter at zero tempera- 
ture and finite chemical potential, the phenomenological 
parameters of quasi-particle model at zero temperature 
and finite chemical potential has some uncertainties. The 
main purpose of this paper is try to constrain the parame- 
ters of the quasi-particle model at finite chemical poten- 
tial based on the most recent astronomical observation 
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This paper is organized as follows. In Section II, we 
introduce the quasi-particle model and its parameters. 
In Section III, the structure of neutron-like star is dis- 
cussed. The relation between the parameters of quasi- 
particle model and mass-radii relation of neutron- like star 
is also illustrated. Through comparison of theoretical 
results with astrophysical observation Q, one can give 
a strong constraint on the parameters of quasi-particle 
model used in the present work. Finally, a summary is 
given in Section IV. 



II. THE QUASI-PARTICLE MODEL AND ITS 
PARAMETERS 

Recently, by means of functional integral formalism a 
model-independent formula for calculating the EOS of 
QCD at finite chemicalpotential and zero temperature 
was proposed in Ref. [34j. By using this formula the 
authors in Ref. (35j obtained an explicit analytic expres- 
sion for the EOS using the quasi-particle model. In this 
paper, we shall use the EOS of this model to study the 
structure of the neutron-like star. The EOS in Ref. [35[ 
is obtained by using path integral method with a effective 
quark propagator. This quark propagator has the form 
of a free quark propagator with a density dependent ef- 
fective mass. 

According to Ref. [34j , the quark number density reads 

p(ji) = —N c Nf J j^HGMiPhih (!) 

where G\p](p) is the full quark propagator at finite chem- 
ical potential fi, N c and Nf denote the number of colors 
and of flavors, respectively, and the trace operation is 
over Dirac indices. According to the above formula, in 
order to calculate the quark number density, one should 
know the full quark propagator at finite /1 in advance. 
Unfortunately, at present no one knows how to calculate 
the full quark propagator from first principles of QCD. 
So, one has to choose some effective model quark propa- 
gator as the approximation of the full quark propagator. 
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In this paper we choose the following model quark prop- 
agator proposed in Ref. |35j 



G VKp) = *7 ■P + m(fi) ) 



(2) 



where p = (p,p4 + i/i) and the Euclidean gamma matrices 
satisfy the algebra {7^,71/} = 25 uu . Here, following Refs. 
[H,|33], we choose the effective quark mass m(/x) and the 
effective coupling constant <? 2 (/i) to be: 



m 2 (/z) = 



mV(m) 

3tt 2 : 



<? 2 (m) = 



I6vr 2 



91n(a(/x + c)) 2 



(3) 



where a, c are phenomenological parameters of quasi- 
particle model. In this work, the current mass of u,d,s 
quarks is taken to be zero. It could satisfy chemical equi- 
librium and neutrality naturally. Then we can use the 
contour integration method to obtain the quark number 
density 



NgNf 

3tt 2 



{v?-m\n)flH{n-m{v)). (4) 



Due to the step function in the right-hand side of Eq. (4) , 
it can be found that the quark number density vanishes 
when \i is below a critical value fiQ (In this work, the 
value of no depends on the parameters a and c) . Namely, 
fx = fj,Q is a singularity which separates two regions with 
different quark number densities. This behavior a gree s 
qualitatively with the general conclusion of Ref. |38| . 
Now let us turn to the calculation of the EOS. According 
to Refs. [H HI, the EOS of QCD at T = reads 



P( M )=P( M )| M=0 



(5) 



where P(fj,)\fj,=0 is vacuum pressure which represents the 
pressure density at /x = 0. In this paper, we take 
-P(/x) Im=o = — where B is a phenomenological param- 
eter in our present work (at present one cannot calculate 
its numerical value reliably from first principles of QCD). 
Here it should be noted that Eq. (5) is a model inde- 
pendent formula. This formula shows that the pressure 
density P{p) at finite \i and zero temperature is totally 
determined by the quark number density p(p) (up to a 
constant term P(p)\ fl= o)- Therefore, the quark number- 
density has the correct behavior required by QCD im- 
plies that the pressure P(n) also has the correct behav- 
ior required by QCD at finite chemical potential and zero 
temperature. 

Substituting Eq. (4) into Eq. (5), we obtain 

p^) = -b + \ f <*y<v - m^'W 2 - m V)) 3/2 , 

n Jo 

(6) 

where we have taken Nf = N c = 3. From Eq. (6) it can 
be seen that the EOS of quasi-particle model depends 
on the model parameters a, a Moreover, the so-called 
vacuum pressure B is also a parameter which affects the 
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- SPressure of QGP and Nuclear matter 
B 1 '"=1 37MeV of QGP 

SPressure of QGP and Nuclear matter 
B 1,4 =140MeV of QGP 

- SPressure of QGP and Nuclear matter 
B 1M =1 50MeV of QGP 
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Fig. 1. The SPressure of quark matter described by 
quasi-particle model (a — 12/GeV,c — —0.05GeV) and 
hadronic matter with different B 



EOS in the study of neutron-like star }4lH43l |. In the fol- 
lowing, we shall discuss the relation between parameters 
a, c and B. 

The choice of B in quasi-particle model is compli- 
cated (depending on the corresponding hadronic EOS 
44j and quark EOS). The different parameters a,c in 
quasi-particle model lead to different ranges of B. In or- 
der to see it clearly, let us to define 8P(n) = Pquarkiv) — 
Phadron (/•*) t° study the phase transition. Here P qua rk (/■*) 
and Phadronin) denote the pressure of quark matter and 
hadronic matter, respectively. P qU ark(p) can be obtained 
using Eq. (6). According to the EOS in Ref. [iij . 
Phadron = when /i < 307MeV. A physical phase 
transition infers that SP(fj-) intersects with once when 
fx > 307 MeV, that is B > Pg« O rfc(0.307GeV%^ , where 
Pquark(p)\ fj.=io stands for the second term in the right- 
hand side of Eq. (6). However, the case is not always 
so simple. We show 5P(/j,) in Fig. I with model param- 
eters a = 12/GeV, c = — 0.05GeV. Just as is shown 
in Fig. I, SP(fi) is not even a monotonously increas- 
ing function. There are two peaks SP(hi) and 8P(p,2)j 
Hi < fJ-2- We have two choices to get a valid result: (I) 
5P(ni) < 0, that is B > P^ark (Mi)U^o - fWfr<m(Mi)> 
i.e., B 1 / 4 = IbOMeV in Fig. 1; (2) SP(fx 2 ) > and 
B > P quark (Q.307GeV), that is P OU a r k(0-307GeV)\^o < 

B < Pquark(H2)\u^0 - Phadron^), i.e., P 1/4 = 137 'M eV 

in Fig. 1. It can be seen that some values of B (i.e., 
P 1 / 4 = 140MeV) would lead to the result which is not 
physical. It may be due to the inconsistency of QGP 
model and constant vacuum pressure. At present we do 
not have a valid method to deal with it, instead, when 
we use QGP model to study the structure of the neutron- 
like star, we should pay attention to the range of B. For 
some sets of model parameters a, c, the second choice 
does not exist because P qU ark{H2)\^o - Phadronini) < 
P q uark{0.307GeV)\^ . 

The mass-radii relation of neutron-like star depends 



3 



-Hybrid B =132.7MeV 
- - Hybrid B"*=142.04MeV 
■Hybrid B'"=149.83MeV 




■ QGP B -132.7MeV 
■■ QGP B ,/4 =142.04MeV 
-QGP B 1M =149.83MeV 
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Fig. 2. The mass-radii relations of hybrid star and pure 
quark star with a — 13/GeV, c = — O.lGeV and different B 

obviously on B. The mass-radii relation of the pure 
quark star and the hybrid star (with a — l3/GeV,c = 
— O.lGeV) with different B are plotted in Fig. 2 . From 
Fig. 2 it can be seen that for pure quark star lager B 
yields smaller mass and radius, while the case for hybrid 
star is more complicated. We will discuss it in Part B of 
Section III. 



III. THE STRUCTURE OF NEUTRON-LIKE 
STAR 

Using the EOS, we can obtain the structure 
of a neutron-like star by integrating the Tolman- 
Oppenheimer-Volkoff equation 



dP{r) G(e + P)(M + ATrr 3 P) 



dr 
dM(r) 
dr 



r(r - 2GM) 



= Aur e. 



The energy density reads 



dP 

e{jM) = -P(n) + n—. 



(7) 
(8) 

(9) 



In this paper, we consider the interior structure of a 
neutron-like star composed, in whole or in part, of quark 
matter. The structure of pure quark star is described to- 
tally by the quasi-particle model. The hybrid star which 
is proposed by the authors of Refs. 0-Q has a quark core 
described by quasi-particle model and a crust of nuclear 
matter described by the APR hadronic EOS 



A. Pure Quark Star 

In this part, we use the recent observation of 1.97 ± 
0.04-solar-mass neutron-like star Q to constrain the pa- 
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Fig. 3. The maximal mass changes with the parameter a for 
different c. 



rameters a, c. In the previous section, we show that the 
larger B yields smaller mass, so we can choose the small- 
est B which leads to the valid result to get the maximal 
mass for every set of a, c. 

The relation between the mass and model parameters 
a, c is shown in Fig. 3. It can be seen clearly that for 
the same c, the mass for different a increases at first, 
then decreases later. This is because at first the smallest 
B w 0, larger a yields larger mass. When B dominates 
the maximal mass, larger a leads to larger B which yields 
smaller mass. This two obvious regions come from the 
different choices of B which we discuss in Section II. For 
different c, larger c yields larger mass at first, then the 
larger c leads to smaller mass. The largest maximal mass 
to a increases with c increasing. 

With the mass related to the the parameters a, c, we 
can constrain the parameters with the recent observation 
of 1.97±0.04-solar-mass neutron-like star Q. We assume 
it is a pure quark star which can be described by quasi- 
particle model. 

The range of a, c is illustrated in Fig. 4. The physical 
meaning of the parameters a, c is not obvious. Here, we 
try to understand it from the point of view of the running 
coupling constant a s . The comparison of the experimen- 
tal data on a s and the results of a s calculated with some 
representative values of a, c are shown in Fig. 5. Here, 
one point should be clarified. In the usual expression of 
the running coupling constant a s (Q 2 ), a s is a function of 
the momentum scale Q. In our work, in order to under- 
stand the physical meaning of a s , the chemical potential 
H is taken as an energy scale similar to the momentum 
scale Q in the expression of a s . From Fig. 5 and Fig. 
6 (see below), it can be seen that the parameters which 
yield 1.97±O.O4M0 are also consistent with the exper- 
imental data on a s (45[. Therefore, we conclude that 
quasi-particle model with proper parameters is suitable 
to describe quark matter. 

If we assume the central density /i of quark star to be 
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Fig. 4. The contour of the parameters a, c constrained by 
recent observation 0] 
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Fig. 6. Mass of different central potential with different 
parameters a, c, where p„ui = 2.8 x f0 14 g / cm 3 . B is chosen 
to be smallest which is physical to obtain the maximal mass. 



experiment data 
- a=5/GeV,c=-0.02GeV ■ 
-a=6/GeV,c=-0.05GeV ■ 
-a=8/GeV,c=-0.1GeV ■ 
a=10/GeV,c=-0.12GeV 
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Fig. 5. Comparison of the running coupling constant a s 

calculated using different values of a, c with the 
experimental data. The experimental data are taken from 
table 1 of Ref. g|. 



0A5GeV, from Fig. 6 we deduce it is reasonable, then the 
range of a s is from 0.659 ~ 1.135. By calculation we can 
also find that the range of [iq is 0.21GeV ~ 0.327GeV. In 
Ref. [HI], based on a universal argument, it is pointed out 
that the existence of some singularity at the point \l = [Ao 
and T = is a robust and model-independent predic- 
tion. The numerical value of critical chemical potential in 
pure QCD (i.e., with electromagnetic interactions being 
switched off) is estimated to be Mcv _ 3Q7 MeV 

(where is the nucleon mass and N c = 3 is the num- 
ber of colors). Here, in order to be self-contained, we 
shall give an explanation for how this result is derived. 
According to the argument of Ref. [11], the critical 
baryon chemical potential ^bo = xnin a (E a /N a ), where 
a stands for a generic quantum state of the system and 
E a and N a is the energy and baryon number of this 
state, respectively. The energy per baryon, E/N, can 



also be written as — (jVmjy — E)/N, where m/v is 
the nucleon mass. Therefore, the state which minimizes 
E/N is that for which the binding energy per nucleon, 
e = (N to 7v — E)/N, is a maximum. Empirically, we 
know that this state is a single iron nucleus at rest with 
N = A = 56 and e w 8 MeV. However, in QCD with- 
out electromagnetism the binding energy per nucleon in- 
creases with N = A. This is the consequence of the 
saturation of nuclear forces and can be seen from the 
Weizsacker formula. Without electromagnetism, only the 
bulk and surface energy terms are significant for large A: 



e(A) 



Attin — m,A 
A 



ai — CI2A 



-1/3 



with at w 16 MeV, a 2 « 18 MeV 46]. As A -> 00, e sat- 
urates at the value a\. This corresponds to the binding 
energy per nucleon in a macroscopically large sample of 
nuclear matter as defined by Fetter and Walecka in (4|| . 
Therefore, in QCD the critical baryon chemical poten- 
tial /iso ~ m N — 16 MeV and the corresponding quark 
chemical potential /j is m "-i6 MeV _ 3Q7 Me y_ The 



range of parameters a, c constrained by the observation 
of 1.97±0.04-solar-mass neutron-like star [3J is consistent 
with it. 



B. Hybrid Star 

The hybrid star was inspired by the idea of the phase 
transition from quark matter to nuclear matter 04a, S3] ■ 
The phase transition in this paper is modeled by a simple 
Maxwell's construction. 

In Fig. 7 we plot the relation of the maximal mass to 
the parameter B. From Fig. 7 it can be seen that in 
the large B region, the maximal mass tends to 2.2Mq. 

This means that for large B, the hybrid star tends to 
a star which is composed totally by hadronic matter. 
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Fig. 7. The maximal mass changes with the parameter B 
with different a, c. 
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The maximal mass for small B depends on the parame- 
ters in quasi-particle model. For some parameters (e.g., 
a = l2/GeV,c = —O.lGeV), the maximal mass in- 
creases with B increasing, while for other parameters 
(e.g., a — 12/GeV, c = — 0.05GeV), the maximal mass 
decreases with B increasing at the beginning, then it in- 
creases with B increasing. This means that for the lat- 
ter parameters, the behavior at small B tends to a pure 
quark star. The discontinuous curves with parameters 
a = 12/GeV, c = -0.05GeV and a = 16/GeV, c = 
— O.lGeV shown in Fig. 7 is caused by the valid ranges 
of B. We can c omp are the result with the same plot pub- 
lished in Ref. [48( which has a quark core described by 
MIT bag model and hadronic crust described by TM1 
RMF EOS and NL3 RMF EOS. However, the running 
coupling constant in the center of the star in Fig. 7 is 
narrow. In Fig. 8, we plot the maximal mass to B with 
different central couplin g co nstant. The plot is similar 
with the result in Ref. [48| , In Ref. [48|], the authors 
stressed that different hadronic model could lead to the 
differences of the result. In our work, we use APR EOS 
to describe the hadronic crust. 

From Fig. 9 and Fig. 11 , a smaller B will narrow the 
region of a-c space. However, for hybrid star we cannot 
give a rigid range of a and c. It can also be figured 
out that smaller a and c will yield larger maximal mass. 
According to Eq. (3), larger running coupling constant 
leads to larger mass, which is consistent with the 1.97 ± 
0.04 solar-mass neutron-like star. 



IV. SUMMARY 

The parameters of quasi-particle model at zero temper- 
ature and finite chemical potential can be constrained by 
recent observation f|. We consider the structure of a 
neutron star which is composed, in whole or in part, of 
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quark matter. In this paper quark matter is described 
by quasi-particle model incorporating a bag constant B. 

We discuss the range of B and the influences of B 
on both quark and hybrid star. First, we assume the 
neutron-like star is a pure quark star. Larger B yields 
smaller mass. We can constrain the parameters a, c by 
recent observation |3j. The range of the running cou- 
pling constant a s with central density /i = 0A5GeV of 
pure quark star is 0.659 ~ 1.135 (this clearly shows that 



even at such a high density the coupling constant is still 
large, so perturbation theory cannot be used). We can 
also find that the range of fi is Q.2\GeV ~ 0.327GW, 
which is consistent with the result given in Ref. [38| . 
Second, we assume the neutron-like star to be a hy- 
brid star which has a quark core and a hadronic crust. 
The parameters of quasi-particle model and the vacuum 
pressure B affect the structure of hybrid star. In the 
large B region, the maximal mass tends to 2.2Mq. 

The maximal mass for small B depends on the parame- 
ters of quasi-particle model. For some parameters (e.g., 
a = 12/GeV,c = — O.lGeV), the maximal mass in- 
creases with B increasing, while for other parameters 
(e.g., a — 12/GeV,c = —0.05GeV), the maximal mass 
decreases with B increasing at the beginning, then it in- 
creases with B increasing. Larger a and c lead to smaller 
mass of hybrid star. We conclude that either a pure 
quark star and a hybrid star with a quark core described 
by quasi-particle model can be consistent with the re- 
cent observation of PSR J1614-2230. The parameters in 
quasi-particle model can be narrowed by the observation. 
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